On the thermal footsteps of Neutralino relic gases 
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Abstract. Current literature suggests that neutralinos are the dominant cold dark matter particle 
species. Assuming the microcanonical definition of entropy, we examine the local entropy per 
particle produced between the "freeze out" era to the present. An "entropy consistency" criterion 
emerges by comparing this entropy with the entropy per particle of actual galactic structures given in 
terms of dynamical halo variables. We apply this criterion to the cases when neutralinos are mosly 
b-inos and mostly higgsinos, in conjunction with the usual "abundance" criterion requiring that 
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present neutralino relic density complies with 0.1 < i2„o < 0.3 for h ~ 0.65. The joint application 
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of both criteria reveals that a better fitting occurs for the b-ino channels, hence the latter seem to 
be favoured over the higgsino channels. The suggested methodology can be applied to test other 
cn : annihilation channels of the neutralino, as well as other particle candidates of thermal gases relics 
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There are strong theoretical arguments favouring lightest supersymmetric particles 
(LSP) as making up the relic gas that forms the halos of actual galactic structures. As- 
suming that R parity is conserved and that the LSP is stable, it might be an ideal candi- 
date for cold dark matter (CDMJl provided it is neutral and has no strong interactions. 
The most favoured scenario d 0. S 0. H 0] considers the LSP to be the lightest neu- 
tralino a mixture of supersymmetric partners of the photon, Z boson and neutral 
Higgs boson |2|]. Since neutralinos must have decoupled once they were non-relativistic, 
it is reasonable to assume that they constituted originaly a Maxwell-Boltzmann (MB) in 
thermal equilibrium with other components of the primordial cosmic plasma. In the 
present cosmic era, such a gas is either virialized in galactic halos, in the process of 
virialization in halos of galactic clusters or still in the linear regime for superclusters and 
structures near the scale of homogeneity [ZLISLbl]- 

The equation of state of a non-relativistic MB neutralino gas is 
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P = m %\ n %\ 1 1 + 2W ' P = ~c~' X = T' (1) 

where m*o and n^o are the neutralino mass and number density. Since we will deal exclu- 
sively with the lightest neutralino, we will ommit henceforth the suscript „o, understand- 
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ing that all usage of the term "neutralino" and all symbols of physical and observational 
variables (i.e. £Iq, m, p, n, etc.) will correspond to this specific particle. As long as the 
neutralino gas is in thermal equilibrium, we have 



m 



1 3 



/27T 



x 3//2 exp (—x) 
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where g = 1 is the degeneracy factor of the neutralino species. The number density n 
satisfies the Boltzmann equation 



n + 3Hn = -(a|v|) n 2 - (n^f 
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where H is the Hubble expansion factor and (o|v|) is the annihilation cross section. Since 
the neutralino is non-relativistic as annilitation reactions "freeze out" and it decouples 
from the radiation dominated cosmic plasma, we can assume for H and (o|v|) the 
following forms 
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where m p = 1.22 x 10 y GeV is Planck's mass, g* = g*(T) is the sum of relativisitic 
degrees of freedom, (v 2 ) is the thermal averaging of the center of mass velocity (roughly 
v 2 oc \/x in non-relativistic conditions) and the constants a and b are determined by the 
parameters characterizing specific annhiliation processes of the neutralino (s-wave or 
p-wave) [2]. The decoupling of the neutralino gas follows from the condition 
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leading to the freeze out temperature Tf. Reasonable approximated solutions of © 
follow by solving for Xf the implicit relation yD 
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where g^ = g*(7f) and cq ~ 1/2 yields the best fit to the numerical solution of © and 
©. From the asymptotic solution of © we obtain the present abundance of the relic 
neutralino gas M\ 
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where So ~ 4000cm 3 is the present radiation entropy density (CMB plus neutrinos), 
p crit = 1.05 x 10- 5 GeVcm- 3 . 



Since neutralino mases are expected to be in the range of tens to hundreds of GeV's 
and typicaly we have Xf- ~ 20 so that 7f < GeV, we can use ~ 106.75 Jj] in equations 
- ©. Equation shows how Xf has a logaritmic dependence on m, while theoretical 
considerations JE H S 0, El 0] related to the minimal supersymetric extensions of 
the Standard Model (MSSM) yield specific forms for a and b that also depend on 
m. Inserting into (EJ)-© the specific forms of a and b for each annihilation channel 
leads to a specific range of m that satisfies the "abbundance" criterion based on current 
observational constraints that require 0. 1 < Clo < 0.3 and h w 0.65 |0]. 

Suitable forms for (o|v|) can be obtained for all types of annihilation reactions J2]. 
If the neutralino is mainly pure bino, it will mostly annihilate into lepton pairs through 
t-channel exchange of right-handed sleptons. In this case the cross section is p-wave 
dominated and can be approximated by <|5]) with 10,0 Ell 
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where m/ is the mass of the right-handed slepton and a 2 = g\/4% ~ 0.01 is the fine struc- 
ture coupling constant for the U(l)y gauge interaction. If the neutralino is Higgsino-like, 
annihilating into W-boson pairs, then the cross section is s-wave dominated and can be 
approximated by © with oLfTolfTTIl 

2m z (2 — m z /m z J z 

where m w = 80.44 GeV is the mass of the W-boson and a 2 = g\l^ — 0.03 is the fine 
structure coupling constant for the SU (2)l gauge interaction. 

In the freeze out era the entropy per particle (in units of the Boltzmann constant k B ) 
for the neutralino gas is given by ILZLISLlSIl 
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where we have assumed that chemical potential is negligible and have used the equation 
of state 0. From Q and (fT2"l) . it is evident that the dependence of Sf on m will will be 
determined by the specific details of the annihilation processes through the forms of a 
and b. In particular, we will use (fTOl) and (fTTT) to compute st from (T7l- (fT2l) . 

After decoupling, particle numbers are conserved and the neutralinos constitute a 
weakly interacting and practicaly collisionless self gravitating gas. This gas initialy 
expands with the cosmic fluid and eventualv undergoes gravitational clustering forming 
stable bound virialized structures |& [E EH E3] • The virialization process involves a 
variety of dissipative effects characterized by collisional and collisionless relaxation 
processes JH, 13, 14]. However, instead of dealing with the details of this complexity, 



we will compare the initial and end states of this gas with the help of simplifying but 
general physical assumptions. 

Consider the microcanonical ensamble definition of entropy per particle for a diluted, 
non-relativistic gas of weakly interacting particles, given in terms of the volume of phase 



space [11311 
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where V and E are local average values of volume and energy associated with a 
macrostate that is sufficiently large as to contain a large number of particles, but suf- 
ficiently small so that macroscopic variables are approximately constant. For a gas 
characterized by non-relativistic velocities v/c <C 1, we have V x l/n <x m/p and 
E oc mv 2 /2 oc m/x. Assuming as the initial and final states, respectively, the decoupling 
(sf, Xf, rif) and the values (5"", x (h> , n <h ') that correspond to a suitable halo structure, the 
change in entropy per particle that follows from (fT3l is 
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where we have used (fT2l to eliminate Sf in terms of Xf. Considering the halo gas as 
a roughly spherical, inhomogeneous and self-gravitating system that is the end result 
of the evolution and gravitational clustering of a density perturbation at the freeze out 
era (the initial state), the microcanonical description is an excellent approximation for 
gas particles near the symmetry center of this system where the density enhancement 
is maximum but spacial gradients are negligible. We will consider then current halo 
macroscopic variables (the end state) as evaluated at the center of the halo: sf, xf, nf. 

Bearing in mind that the density perturbations at the freeze ot era were very small 
(8n f /n f < 1(T 4 , 0, 0, H), the density nf is practicaly homogeneous and so we can 
estimate it from the conservation of particle numbers: = hq (1 + Zf) 3 , and of photon 
entropy: g^fSf = g*oSo (1 + Zf) 3 , valid from the freeze out era to the present for the 
unperturbed homogeneous background. Eliminating (1 +Zf) 3 from these conservation 
laws yields 
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where g* = g 
po/pf and p 
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2.7 K. Since for present day conditions no/n 
p cr j t Q.oh 2 , we collect the results from ( [131 and write (fT4l) as 
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Therefore, given m and a specific form of (o|v|) associated with a and b, the entropy 
per particle of the neutralino halo gas depends on the initial state given by Xf in © and 
(fT2l . on observable cosmological parameters £Iq, h and on state variables associated to 
the halo structure. 



If the neutralino gas in present halo structures strictly satisfies MB statistics, the 
entropy per particle, s®, in terms of = mnf and = mc 2 / (k B T^), follows from 



the well known Sackur-Tetrode entropy formula [ 15 
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Such a MB gas in equilbrium is equivalent to an isothermal halo if we identify ifl^l 



1 _ k B T' h> _ 2 
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where is the velocity dispersion (a constant for isothermal halos). 

However, an exacly isothermal halo is not a realistic model, since its total mass di- 
verges and it allows for infinite particle velocities (theoreticaly accessible in the ve- 
locity range of the MB distribution). More realistic halo models follow from "energy 
truncated" (ET) distribution functions 
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that assume a maximal "cut 

off" velocity (an escape velocity). Therefore, we can provide a convenient estimate of 
the halo entropy, s®, from the microcanonical entropy definition (IT~5T> in terms of phase 
space volume, but restricting this volume to the actual range of velocities (i.e. momenta) 
accessible to the central particles, that is up to a maximal escape velocity v e (0). From 
theoretical studies of dy namical and thermodynamical stability associated with ET dis- 
tribution functions II 1 71 1 1 8l 1 1 9l.l22l.l2 lL l23l 120 1 and from observational data for elliptic 
and LSB galaxies and clusters ll24i l25l ud l27l l28h , it is reasonable to assume 



v 2 (0) = 2|<I>(0)| ~ aa 2 (0), 12<a<18, 



(20) 



where 4>(r) is the newtonian gravitational potential. We have then 
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where we used4!' = c 2 /o 2 (0) as in ( fT9l) . As expected, the scalings of (OTT) are identical 
to those of (fTSl . Similar entropy expressions for elliptic galaxies have been examined in 
& 

Comparison between obtained from (OTT) and from (fTTT) leads to an "entropy con- 
sistency" criterion. Since d2TT) scales with In m 4 , while (fTTT) does so approximately with 
In m 3 , we have a weak logarithmic dependence of on m. Therefore, the fulfilment of 
the "entropy consistency" criterion identifies a specific mass range for each dark mat- 
ter particle. This allows us to discriminate, in favour or against, suggested dark matter 
particle candidates and/or annihilation channels by verifying if the standard abundance 
criterion ® is simultaneously satisfied for this range of masses. It is interesting to no- 
tice that both equations, (IPTl and (ED, display an identical functional dependence of sj? 
on the present day halo parameters, p£' and jc£\ This implies that a given dark matter 
particle candidate, characterized by m and by specific annihilation channels given by Xf 
through 0, will pass or fail to pass this consistency test independently of the details one 
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FIGURE 1. Figures (a) and (b) respectively correspond to the higgsino and b-ino channels. The figures 
display sf as a function of log 10 m, obtained from d20i — d2 II (gray strip), from Jl 8I > (crosses) and from d!7i 
for h = 0.65 and Q.q = 0.2 (thick curve), together with its uncertainty strip Q.q = 0.2 ±0.1. The vertical 
strip marks the range of values of m that follow from (|8|l-(|9|l for the same values of ilo and h. It is evident 
that only the b-ino channels allow for a simultaneous fitting of both the abundance and the entropy criteria. 



assumes regarding the present day dark halo structure. However, the actual values of s% ] 
for a given halo structure, whether obtained from (OTT) or from (flTt . do depend on the 
precise values of and x^. 

We will now evaluate d2lT) and (fTTt for the two cases of neutralino channels: the 
b-ino and higgsino, including numerical estimates for X th ' and p' h| that correspond to 
central regions of actual halo structures. Considering terminal velocties in rotation 
curves we have v^ erm ~ 2o^ (0), so that xj? ~ 2(c / 'vterm) 2 > while recent data from LSB 
galaxies and clusters S HH HO, EH H2] suggest the range of values 0.01M Q /pc 3 < 
p* 1 < lM /pc 3 . Hence, we will use in the comparison of (fTTt and (|2TT) the following 
numerical values: p£' = 0.01 M /pc 3 = 0.416 GeV/cm 3 andx™ — 2x 10 6 typical values 
for a large elliptical or spiral galaxy with vterm — 300 km/sec ll3QL l3ll 13211 . Figure la 
displays the as a function of log 10 m, for the halo structure described above, for the 
case of a neutralino that is mostly higgsino. The shaded region marks given by (OTT) 
for the range of values of a, while the vertical lines correspond to the range of masses 
selected by the abundance criterion © for Q.q = 0.1, 0.2, 0.3. The solid curves are 
given by ( fTTt for the same values of ^o, intersecting the shaded region associated with 
(|2TT) at some range of masses. However, the ranges of coincidence of a fixed (fTTt curve 
with the shaded region (I2TT) occurs at masses which correspond to values of that are 



different from those used in (TP71) . that is, the vertical lines and solid curves with same 
intersect out of the shaded region. Hence, this annihilation channel does not seem to 
be favoured. 

Figure lb depicts the same variables as figure la, for the same halo structure, but 
for the case of a neutralino that is mostly b-ino. In this case, both the abundance 
and the entropy criterion yield consistent mass ranges, which allows us to favour this 
annihilation channel as a plausible dark matter candidate, with m lying in the narow 
ranges given by this figure for any chosen value of ^o. As noted above, the results of 
figures la and lb are totaly insensitive to the values of halo variables, x'*' and p®>, used 
in evaluating (|2lT) and (IPTl) . Different values of these variables (say, for a different halo 
structure) would only result in a relabeling of the values of along the vertical axis of 
the figures. 

We have presented a robust consistency criterion that can be verified for any annihi- 
lation channel of a given dark matter candidate proposed as the constituent particle of 
the present galactic dark matter halos. Since we require that of present dark matter 
haloes must match derived from the microcanonical definition and from freeze out 
conditions for the candidate particle, the criterion is of a very general applicability, as it 
is largely insensitive to the details of the structure formation scenario assumed. Further, 
the details of the present day halo structure enter only through an integral feature of the 
dark halos, the central escape velocity, thus our results are also insensitive to the fine 
details concerning the central density and the various models describing the structure 
of dark matter halos. A crucial feature of this criterion is its direct dependence on the 
physical details (i.e. annihilation channels and mass) of any particle candidate. 

We have examined the specific case of the lightest neutralino for the mostly b-ino 
and mostly higgsino channels. The joint application of the "entropy consistency" and 
the usual abbundance criteria clearly shows that the b-ino channel is favoured over the 
higgsino. This result can be helpful in enhancing the study of the parameter space of 
annihilation channels of LSP's in MSSM models, as the latter only use equations © 
and (El)-© in order to find out which parameters vield relic gas abbundances that are 
compatible with observational constraints lllLl2ll2iyiblld]- However, equations © and 
®-© by themselves are insufficient to discriminate between annihilation channels. 
A more efficient study of the parameter space of MSSM can be achieved by the joint 
usage of the two criteria, for example, by considering more general cross section terms 
(see for example [2]) than the simplified approximated forms (fTOI) and (ITTb . This work 
is currently in progress. 
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